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Introduction: 


We  consider  a  body  to  have  an  initial  stress  field  ^  When  the  body  is  loaded  beyond 
the  elastic  range, Aplastic  strain  e^  occurs  in  some  part  of  the  body.  If  the  body  is  then 
unloaded,  the  plastic  strain  remains  and  causes  a  residual  stress  field. T*,  Now  consider  the 
load  applied  again.  This  applied  load  causes  a  stress  t^.  There  is  no  increase  of  plastic 
strain  during  this  process  of  unloading  and  reloading.  ^Now  the  stress  field  is  the  sum  of  the 
initial,  the  residual  and  the  applied  stress  / 


T-  ■  =  r.I  4.  r-8  +  rA 

*”lj  Lij  T  lij  T  ‘•Jj 


(1) 


As  shown  by  Lin,  1968,  using  this  analogy  between  plastic  strain  and  applied  force,  we  can 
find  in  terms  of  e^  in  the  body.  As  the  load  is  increased,  t*  increases.  Then  from  the 
incremental  stress  strain  relation  of  the  material,  the  incremental  plastic  strains  in  different 
parts  of  the  body  are  found  (Lin,  1968).  In  this  inelastic  analysis,  the  relation  between  the 
plastic  strain  and  the  residual  stress  is  needed.  Plastic  strain  is  one  type  of  inelastic  strain. 
Hence  the  reciprocal  relation  of  inelastic  strain  and  residual  stress  is  very  useful  in  numerical 
inelastic  analyses.  -  '  .  Ti'  ^ 
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Derivation  of  the  Reciprocal  Relation  for  Inelastic  Bodies: 


Referring  to  a  set  of  rectangular  coordinates  ( x t ,  x:.  x?),  the  equilibrium  condition  is 
ctven  as 


Tijj  +  F,  =  0 


in  the  interior  of  the  volume  V',  and 


(2) 


"'■ij  Vj 


(3) 


on  the  boundary  F  where  is  the  stress  component.  F,  the  body  force  per  unit  volume  along 
x, -direction.  S,  is  the  traction  per  unit  area  along  x,-axis  and  Vj  is  the  direction  cosine  of  the 
extenor  normal  of  the  boundary  with  Xj-axis.  The  repetition  of  the  subscript  denotes  summa¬ 
tion  from  one  to  three,  and  the  subscript  after  comma  denotes  differentiation  with  respect  to 
that  axis.  Multiplying  (2)  and  (3)  by  a  virtual  displacement  Uj  and  integrating  yields 

J  Ft  Uj"  civ  +  j  S,  Uj'  dr  =  -  J  Tijj  Ui'  dv  +  j  T,j  Vj  u  ‘  dr 

=  1  Xj  »i"j  *  =  J  Xj  =.j  dv 

The  stress-strain  relation  of  an  elastic  body  may  be  written  as 


^ij  —  cijkl  eld 

where  c.jkJ  is  the  elastic  constant. 

j  x.j  ej  dv  =  J  cljkl  ew  e,j  dv  =  j  e^  dv 
From  (4)  and  (6) 

J  F,  u*  dv  +  j  S(  uf  dr  =  J  Fj*  u,  dr  +  j  S'  u,  dr 
This  is  Betty's  reciprocal  theorem  for  elastic  bodies. 
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(5) 


i 


i . 


vs 


I  or 


Now  consider  the  presence  of  inelastic  strain  etJ  which  may  be  a  combination  of  plas¬ 
tic  strain,  creep  strain  and  thermal  strain.  The  total  strain  e,j  is  composed  of  the  elastic  strain 
e,j  and  inelastic  strain  e,. 


Tj  —  cijkl  (^kl  ekl) 

Substituting  this  into  (2)  and  (3)  yields 

Cijkl  (®ldj  —  '-klj)  +  Fj  —  0 


(8) 

(9) 


Sj  —  c tjkJ  (ekl  ekll  vj 

It  is  seen  that-cijlk.  et;  j  and  cljiti  Vj  have  the  same  effect  on  the  strain  distribution  Cjj  as  F, 
and  S,  respectively  and  hence  they  are  denoted  respectively  as  Fj  and  Sj  (Lin,  1968). 


Now  consider  an  inelastic  strain  e^1  in  region  I  giving  F/  and  Sj.  The  displacement, 
strain  and  stress  in  another  region  II  caused  by  e,^  are  denoted  by  u/1,1  e,-1,1  and  x,-1,1  respec¬ 
tively  as  shown  in  Fig.  1. 
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Applying  Betty’s  Reciprocal  Theorem  (7)  to  Figs.  1  and  2  gives 

Jv,  Ll  u,U1  dv  +  /p,  §1  uiUI  dr 

=  jvn  L11  u/W  dv  +  Jrfl  S/1  u*1'1  dT 

Replacing  F,  by  -  cijkJ  e^j  and  S,  by  cljkl  e^  Vj  yields 

-|v,  cijU  e^j  u,UI  dv  +  j  cm  c3  vj  uiUI  dr 

r* 

=  -  Jvn  Cijki  ew”  U,11-1  dv  +  cm  e^11  Vj  u,1-11  d T 
Using  divergence  theorem  and  noting  c^  =  Cyy  Eq  (11)  becomes 
J,,  CijU  Ca  e,j'"  dv  =  /_,,  CijU  ej1  e'J1  dv 

L^'eu'dv^Lxa-'eJ'dv 


For  the  case  v1  =  v11  and  eu‘  *  e^1 ,  the  relation  (12)  still  holds. 

if  is  the  stress  in  region  I  causes  by  inelastic  strain  e^1  in  region  II  and  T,j  ’  is  the 
stress  in  region  II  caused  by  inelastic  strain  e^1  in  region  I.  If  regions  I  and  II  have  no  over¬ 
lap,  the  stress  x^11  and  T^J'1  are  residual  stresses.  If  there  is  overlap  as  shown  in  Fig.  3.  The 
residual  stress  in  the  overlapped  region,  denoted  hy  v1  n  v11,  is  then 

-I.II_-I.II_..  "ii 
'•Id,  '•Id  ''klmn  cmn 


^ki,  ^kl  cklmn  emn 


This  gives 


i  ^  dv  -  f ,  _  .  .  e 3  dv  +  fy,  ^  v«  Til"  dv 
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Fig.  3.  v1  and  v11  have  an  overlap  v1  n  v2 


Application  to  Fatigue  Band  Calculations: 

A  face-centered-cubic  crystal  has  lour  slip  planes,  on  each  of  which  there  are  three  slip 
directions  giving  twelve  slip  systems.  Slip  in  a  slip  system  depends  on  the  resolved  shear 
stress  and  is  independent  of  the  normal  stress  on  the  sliding  plane  (Taylor,  1938).  Let  us  con¬ 
sider  a  most  favorably  oriented  crystal  in  a  F.C.C.  polycrystal.  When  this  polycrystal  is 
loaded,  slip  may  occur  in  some  region  of  this  crystal.  After  unloading,  this  slip  remains  and 
causes  a  residual  stress.  Imagine  that  the  polycrystal  is  reloaded.  This  resolved  shear  stress 
in  this  crystal  can  be  expressed  as  the  sum  of  the  initial  resolved  shear  stress  x1,  the  resolved 
shear  stress  due  to  the  applied  loading  xA  and  the  residual  shear  stress  x^.  If  this  resolved 
shear  stress  reaches  the  critical  shear  stress  xc,  i.e. 

t  =  ta  +  TA  +  tR  =  Tc  ([7) 

in  some  region,  this  region  will  start  or  continue  to  slide. 

We  consider  three  thin  slices  in  a  most  favorably  oriented  crystal  at  the  free  surface  of 
the  polycrystal  as  shown  in  Fig.4  (Lin,  1974).  An  initial  tensile  strain  e,^  is  assumed  to  exist 
in  the  slice  R.  This  e^a  causes  an  initial  compressive  stress  in  R,  a  positive  resolved  shear 
stress  x^p  in  P  and  a  negative  x^  in  Q  (Lin  and  Lin,  1988).  A  cyclic  loading  of  tension  and 
compression  on  the  polycrystal  causes  alternate  sliding  in  P  and  Q  causing  the  growth  of  an 
extrusion  or  an  intrusion  (Lin,  1977,  Lin  and  Ito,  1969).  An  extrusion  causes  a  tensile  strain 
and  hence  a  tensile  stress  along  extrusion  direction.  This  can  cause  the  resolved  shear  stress 
in  a  second  slip  system  to  slide  in  R.  This  secondary  slip  system  can  slide  concurrently  with 
this  primary  slip  system  in  P  and  Q.  Let  a(3  denote  the  primary  slip  system  and  ^T),  the 
second  slip  system.  The  plastic  strains  caused  by  slip  in  these  two  systems  are  denoted  by 
e£p  and  ef^  respectively.  The  resolved  shear  stress  x£p(x)  in  a(3-system  at  point  x  caused  by 
a  unit  e£p  at  x  in  a  differential  volume  dv  is  expressed  as  G(x,  aP;  x,ap)e„(j  dv.  If  e^  also 
occurs,  we  can  write 


7 


Nil 


Most  Favorably 
Oriented  Crvstai 


rree  Surface 


PRQ 


*a  (Along  Slip  Direction) 


Grain  Boundary 


mm 


Fig.4  Most  Favorable  Oriented 
Crystal  at  Free  Surface 
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tJjCx)  =  j  G(x.  a(3;  x.  ap)  e^(x)dv  +  J  G(x,  aP;  x.  ^rj)  e^Cx)  dv 


The  two  integrals  are  integrated  over  the  regions  with  plastic  strain.  Similarly  w 


(19) 


e  can  write 


Tc^(x)  =  j  G(x,  qrj;  x.  aP)  e^(x)  dv  +  J  G(x,  qt);  x,  qn)  e^fx)  dv 


(20) 


In  numerical  calculation,  the  thin  slices  P .  Q,  R  are  divided  into  grids.  In  each  grid,  the 
e^jj  or  e^,  is  generally  assumed  to  be  constant  or  linear.  The  residual  resolved  shear  stress 
due  to  plastic  strains  in  the  nth  grid  is  written  as 

t<5j(x)  =  G(x,  ap;  n,  a{3)  e^ov0  +  G(x,  ap;  n,  4r|)  ef,  vn 


where 


J  G(x.  ap;  x,  a(i)  e£p  dv 

G(x,  ap;  n,  ap)  =  — - - - - - 

J  eapOD  dv 

Va 


(22) 


and  e£(in  is  the  average  e<^  in  the  nth  grid.  Now  we  consider  the  residual  stress  in  the  mth 
grid  corresponding  to  an  assumed  plastic  strain  distribution  in  the  grid.  As  mentioned  before, 
this  distribution  is  generally  assumed  to  be  constant  or  linear.  Multiplying  (21)  by  e„p  dv  in 
this  grid  and  integrate 


\  T5a(x)e£fj(x)dv  =  e£0ovn  J  G(x,  ap;  n.  ap)  e£3(x)dv 


+  efi„  vn  /  G(x,  aP;  n,  aP)  e^0(x)  dv 

vm 

=  ea0„  vn  G(m,  n,  ap)  e£jjB  vn  +  e^  vra  G(m,  aP;  n,  ^t|)  e^m  vm 


(23) 


where 


G(m.  a|3;  n.  a(3)  = 


(24) 


Jv  G(x,  ap;  n.  ap)  e^(x)dv 
L  e«0(x)dv 

v  m 

Writing  the  resid-  al  resolved  shear  stress  in  the  mth  grid  as  the  left  hand  side  of  (23)  divided 
by  eniimvr.  ~  j  ea|j(x)  dv,  we  have 

vai 

T&a  =  G(m’  aP;  n-  aP)  ea3o  vn  +  G(m,  ap;  n,  %r\)cL  vn 

(25) 

Letting  the  region  I  be  the  nth  grid  and  region  II,  the  mth  grid,  the  left  hand  side  of  (23) 
becomes  the  left  hand  side  of  (12).  Hence  the  reciprocal  relation  holds.  Similarly  we  can 
write 

t£u  =  G(m-  n,  ap)  e^8vn  +  G(m.  qrj;  n,  %r\)  e^o  vn 

(26) 

G(m.  aP;  n.aP)  and  similar  terms  are  the  influence  coefficients  for  residual  stresses  caused  by 
inelastic  strains  in  different  grids.  Eq.  (16)  gives  the  reciprocal  relation  of  these  influence 
coefficients.  If  the  plastic  strain  is  taken  to  be  constant  in  each  grid,  Eq.  (22)  reduces  to 

G(x,  ap;  n,  ap)  =  —  |v  G(x,  aP;  x,  ap)  dv 

vn  3 

and  Fig.  (24)  reduces  to 


G(m.  aP;  n.  aP)  =  —  [  G(x,  aP;  n,  aP)dv 

V  J 


If  there  are  M  grids  with  incremental  Ae£,  and  N-gnds  with  Ae^.  The  residual 
stresses  are  written  as 


AT&ra  -  £  G(m.  aP;  n.  §r\)  Ae^v,,  +  £  G^m’  aP;  n-  aP)  Ae&n  vn 
n=l  n=l 
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(28) 


i>l  IN 

=  Z  G  (m-  n-  5n)  AeJ,  Vn  +  XG(m,  ^T  n,  a(3)  Ae£Pnvn 

n=l  n=l 

In  finding  the  incremental  plastic  strains  in  different  grids,  solution  of  a  system  of  simultane¬ 
ous  equations  is  required.  Due  to  the  reciprocal  relation,  the  matrix  of  these  equations 
become  symmetrical.  This  helps  greatly  the  numerical  solution. 
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